Cosmological constraints on the Hu-Sawicki modified gravity scenario 
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In this paper we place new constraints on a f{R) modified gravity model recently proposed by 
Hu and Sawicki. After checking that the Hu and Sawicki model produces a viable cosmology, i.e. a 
matter dominated epoch followed by a late-time acceleration, we constrain some of its parameters by 
using recent observations from the UNION compilation of luminosity distances of Supernovae type 
la, including complementary information from Baryonic Acoustic Oscillations, Hubble expansion, 
and age data. We found that the data considered is unable to place significant constraints on the 
model parameters and we discuss the impact of a different assumption of the background model in 
C ' cosmic parameters inference. 
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I. INTRODUCTION 



The recent cosmological data from Cosmic Microwave Background Anisotropics, galaxy surveys and luminosity 
^ ^1 distance of type la supernovae are all providing supporting evidence for a dark energy component, responsable for 
I ■ more than 70% of the total energy budget in our universe (see e.g. [l|). Several candidates have been suggested 
O ■ for explaining this component, as, for example, minimally coupled scalar fields (see e.q. 0] and references therein). 
> I However, it may also be that the cosmological evidence for acceleration comes from the wrong assumption of general 
relativity, i.e. that no dark component is present but actually a modification to gravity is at work. In this respect, 
f{R) theories seem to provide a quite large number of viable models (for a recent review see Q and Q). A particular 
f(R) model that evades solar system test has been proposed by Hu and Sawicki ([5.], HS hereafter). The model has a 
modified Einstein-Hilbert action: 



S = J d + L„,\ (1) 

Lm is the matter lagrangian, — SttG and 

/(i?) = -m \'^^l (2) 
1 + C2(;^) 



?H with m = k p/3 and ci, C2 and n as free parameters. 



As shown in [5] this model is able to reproduce the late time accelerated universe but with distinctive deviations 
from a cosmological constant. In this paper we investigate the cosmological viability of the HS model in more detail. 
After a brief description of the model, in the next section we will show that the HS model satisfies indeed the general 
conditions presented by [6'| as a viable f{R) model. In Sec. Ill we compare the HS model with current data from SN-Ia 
luminosity distances from the UNION catalog (^Sj); Baryonic Acoustic Oscillation data ([HI) and age constraints from 
the analysis of Simon, Verde and Jimenez ([0]). We show that the current data is fully compatible with the HS model 
and that, unless a prior on the matter density is used, the parameters of the model are unconstrained. In particular, 
we analyze the impact of the HS model in the determination of the current matter density. As we will show, assuming 
the HS model instead of the standard cosmological model, could relax the constraints on the effective matter density. 
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A future incompatibility between the values of the matter density determined from different datasets and under 
the assumption of the standard ACDM model could therefore provide an hint for a modified gravity scenario. 

II. THE HU-SAWICKY MODEL 

Let us briefly review in this section the basic equations and results of the HS model. Varying the action in Eq. 1 
with respect to the metric g'^'^ one obtains the modified Einstein equations: 

Gm. + IrR^-^ - (| - - v^v,/fl = er^, (3) 

where fn = df /dR and fj^ji = d^f/dR^ and assuming a flat FRW metric, the modified Friedmann equation: 

- fniHH' + H') + ^ + H'JrrR' = ^ (4) 
6 6 

with ' = d/dlna and p the matter density at present time. 

Defining the new variables uh ~ {H"^ /m?) — and = {R/m?) — 3a~^ the Friedmann equations can be expanded 
in a system of two ordinary differential equations: 

v'r = 9«"' - I -3 2, [VH - fni^ - - ^) + (6) 

In order to compare the HS model with the cosmological constraints usually derived under the assumption of dark 
energy, it is useful to introduce an effective dark energy component with present energy density ri^r = 1 — and 
equation of state Wx, where D,„i is the effective matter energy density at present time. 

Of course, in reality no dark energy component is present, the only component present is matter and modified gravity 
gives the acceleration. Considering the Friedmann equation: 

|! = % + f,,e/.^-^I^ (7) 
the effective equation of state parameter for the dark energy component is given by 

w.^-i-l-y^ (8) 

The free parameters Ci and C2 that appear in Eq. 2 can be expressed in function of the effective density parameters 
by: 

Cl , ^ 

— w 6^ (9) 

C2 

C2 n fl^ 

These relations show that the free parameters of the model are flm, n, and /r^. The latter is constrained to 
I/-R0I ^ 0.1 by solar system tests Q and we will not investigate larger values in the next sections. 
Solving the differential equations system for different values of n, /r^ and f2m, it is possible to obtain various 
evolution trends for the equation of state parameter Wx- In FigUand Fig[2]we plot the behavior of Wx in function of 
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FIG. 1: The equation of state parameter for n — 1 and fn^ equal to —0.1 (solid line), —0.03 (dashed line) e —0.01 (dotted 
line), tlm is set to 0.24 

the redshift z for different values of fn^ and n. 

As we can see in Figll]and Fig[2]the equation of state parameter Wx follows a peculiar behavior in function of the 
redshift. At the present time {z = 0) Wx has always a value higher than the one predicted by the ACDM model 
(w = —1) and, moving towards higher redshifts, it decreases crossing into the phantom region, i.e., assuming values 
lower than —1. For even higher redshifts, Wx moves asymptotically towards —1. 

The same behavior is shown for any value of n and //j^, moreover decreasing the absolute value of f^^ brings Wx 
closer to —1, while increasing n shifts the phantom crossing at lower redshift. 




Z 



FIG. 2: The equation of state parameter Wx for n = 2 and fn^ equal to —0.1 (solid line), —0.03 (dashed line) e —0.01 (dotted 
line). Clm is set to 0.24 

The modification of the Einstein- Hilbert lagrangian brings to a new equation for the expansion of the Universe. 
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The predicted expansion must be consistent with standard cosmological results, i.e., should produce an accelerated 
era after radiation and matter dominance. 

Modified gravity models consistent with current observations, for example, should not change the scale factor evolution 
during the matter era. Hence, it is possible to derive general conditions for the cosmological viability of f{R) theories. 
Introducing the parameters 



m{r) 



RIrr. 



R{l + fR) 



l + fn R + f 

it is possible to show that for f{R) theories the following conditions apply: 
• The model has a standard matter era with a following accelerated phase if 



(11) 



m(r) « and m'{r) > —1 with r w —1 



(12) 



• The accelerated phase goes asymptotically towards the one produced by a dark energy with equation of state 
parameter w = — 1, if 



< m{r) < 1 for r — —2 
• The expansion is not of the phantom type {w < —1) if 



(13) 



m(r) = — 1 — r 



(14) 




FIG. 3: Plot of the two solution for m(r) obtained setting n — 1, Clm = 0.3 and flA = 0.7. The solid hne corresponds to the 
viable solution while the dashed line lies outside the viability region. The red line shows the solution m(r) = — r — 1. 



It is possible to calculate m(r) for the Hu and Sawicki model and to show the cosmological viability of this model. 
In Figl3]we show that, for example, setting n = 1, (im = 0.3 and ^Ia — 0.7, one obtains two solutions for m(r), one 
living outside the viability region and the other corresponding to an acceptable expansion. 
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III. CONSTRAINTS ON THE HS MODEL 



A. Method 



In order to constrain the free parameters of the Hu and Sawicki model {flm, n and /i?o)i we predicted the expected 
theoretical values for a set of observables. 

As now common in the literature, we considered the luminosity distance, defined by: 



dda) = - f 

^ J a 



da ^ f n K\ 



and the Hubble parameter: 



Hia)^ ^nr^HliijH + a-^) (16) 
Moreover, we also considered the quantity: 



A = 



(17) 



where — 0.35, r(z*) — J^* dz/e{z) and e(z) = H{z)/Hq. The value of this parameter can be obtained from 
observations of Baryon Acoustic Oscillations (BAO) 7]. Hence, we have another way to compare model prediction 
with data. 



We used the superovae data from Kowalski et al. J8| to obtain the observational trend of dL{z) and we considered 
H{z) values obtained by Simon, Verde and Jimenez [9| and a prio r on the Hubble parameter Hq = 0.72 ±0.08 derived 
from measurements from the Hubble Space Telescope (HST, flO]), Finally, we used the value of A from Eisenstein et 

al. [m. 

We compute a variable for each observational quantity and then combine the results in a single variable = 
XsN + Xbao + Xh + X%ST- Once the theoretical evolution of the observational quantities is defined, we can define a 
likelihood function as a function of n and fn^ as 



L 



(18) 



where Xmin ^^'^ minimum value in the considered range of n and fug . 



B. Results 



Combining the results obtained from the comparison between the experimental data for H{z), A and dL{z) and 
their theoretical values, we can constrain the free parameters n and //Jq for different values of Vtm and Vlx- 

Setting — 0.2 and ^Ix = 0.8 it is possible to find an upper limit on n and on n < 1.6 and fn^ < —0.03 at 
2 tr, while, performing the same analysis with different values of f2„i and r^A, we can see that raising i7„j brings to 
more loosely constrained parameters. 

We can note anyway that for higher values of flm, higher n are preferred, while smaller values of n are more in 
agreement with data for smaller flm- 

As we can see from Figure |4l for fi™ — 0.3 both parameters are almost totally unconstrained; this points out the need 
of an independent measurement of the effective matter content (flm) in order to obtain some constraints on n and /i?,,. 
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FIG. 4: 68%, 95% and 99% confidence levels in the n-fRg plane in function of different values of 



It is interesting to compare the best fit values of obtained in the modified gravity framework with the of a 
cosmological constant model. In order to quantify the goodness-of-fit of the two models we use the Akaike information 
criterion (AIC) and the Bayesian information criterion (BIC) [l^, defined as 

AIC = -2\nL + 2k 

BIC = -2\nL + k\nN 
where L is the maximum likelihood, k the number of parameters and N the number of points. 
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FIG. 5: AIC (left panel) and BIC (right panel) tests in function of Sim for the standard case of a cosmological constant (dashed 
line) and for the HS model (solid line). 



In Fig. Owe plot the best fit values of the AIC and BIC tests in function of flm for the standard model based 
on a cosmological constant and the HS model respectively. As we can see, while the cosmological constant gives 
slightly better values for the overall best fit, when larger or smaller values of flm are considered the AIC and BIC 
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tests provide definitly better values for the HS model. In few words, there is a weaker dependence of the observables 
considered from f7™ in the case of HS scenario. 



Sn+BAO Sn+H+H, 




0.20 0.25 0.30 0.35 



FIG. 6: The likelihood function L{Qrn) for the Hu and Sawicki model (thicker lines) and AC DM model (fainter lines). 

It is therefore important to quantify the impact of a different choice of the theoretical background model on the 
derived constraints on Qm- In Fig- [6]we compare the constraints on the D,„i parameter derived under the assumption 
of the HS scenario with the similar constraints but assuming general relativity and dark energy. As we can see 
the dm parameter is less constrained respect to the ACDM scenario. This is certainly due to the larger amount of 
parameters present in the HS model. In the future, with the increasing experimental accuracy, if a discrepancy between 
independent constraints on the matter density will be found then a modified gravity scenario could be suggested as 
possible explanation. This result anyway shows that one should also be extremely careful in considering the current 
cosmological constraints because of their model dependence. 

IV. CONCLUSIONS 

In this paper we have compared a modified gravity scenario, the HS model, with several current cosmological 
datasets. We have found that the model is in excellent agreement with recent SN-Ia, BAO and H{z) data. Moreover, 
the parameters of the model are substantially unconstrained by the data considered. This has important effect on 
the current constraints on some parameters as the matter density. We have shown that the assumption of the HS 
model enlarges the current constraints on this parameter by ~ 30%. If a discrepancy between two experimental 
determinations of the matter density will be found in the framework of general relativity, then a possible solution 
could be the introduction of a modified gravity scenario. It will be duty of future experiments to scrutinize this 
interesting possibility. 
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